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Preface 

First-round,  fire-for-effect  artillery  fire  requires  accurate  knowledge  of  the 
meteorological  parameters  affecting  the  trajectory;  namely,  vertical  profiles  of 
temperature,  wind,  and  density.  Traditionally,  these  profiles  have  been  mea¬ 
sured  using  the  balloon-borne  radiosonde.  The  radiosonde  is  accurate  and  has 
been  extremely  well  tested  over  the  last  60  years,  but  it  also  has  serious  disad¬ 
vantages  due  to  its  lack  of  timeliness  and  long  logistical  “tail”  (because  it 
requires  a  helium  or  hydrogen  supply  for  its  balloons). 

The  Mobile  Profiling  System  (MPS)  described  herein  can  overcome  these  tempo¬ 
ral  and  logistical  problems  by  automated,  near  real-time  remote  sensing  methods. 
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Executive  Summary 


A  mobile  atmospheric  profiling  system  has  been  developed  by  the  Battlefield 
Environment  Division  (BED)  of  the  U.S.  Army  Research  Laboratory  (ARL)  and 
the  Environmental  Technology  Laboratory  (ETL)  of  the  Environmental 
Research  Laboratories,  National  Oceanic  and  Atmospheric  Administration 
(NOAA).  This  system,  which  is  capable  of  probing  the  atmosphere  fi’om  the  sur¬ 
face  to  altitudes  over  30  km,  combines  ground-based  instruments,  including  a 
five-beam,  924  MHz  radar  wind  profiler;  a  Radio  Acoustic  Sounding  System 
(RASS);  and  two  passive  microwave  sounders  (measuring  virtual  temperature 
and  moisture  parameters,  respectively)  with  a  receiver  and  processor  for  meteo¬ 
rological  satellite  images  and  soundings. 

This  report  describes  the  method  for  merging  data  fi-om  satellite  and  ground- 
based  remote  sensing  systems  and  presents  results  of  early  tests  of  the  Mobile 
Profiler  System  (MPS),  including  information  on  accuracies  achieved  in  its 
measurements.  Examples  of  some  of  the  many  data  formats  produced  by  the 
MPS  are  also  shown. 

Software  in  the  MPS  produces  profiles  fi'om  the  surface  to  the  highest  satellite 
sounding  level  by  combining  profiles  generated  fi'om  the  suite  of  ground-based 
sensors  with  those  fi’om  a  meteorological  satellite.  The  algorithms  generate 
soundings  of  temperature,  humidity,  wind  velocity,  and  other  meteorological 
variables. 

The  MPS  combines  the  capabilities  of  several  types  of  sensing  systems  to  pro¬ 
vide  atmospheric  soundings  with  a  rapid  refresh  rate  that  greatly  reduces  errors 
caused  by  time  staleness.  The  ability  of  the  ground-based  instruments  in  the 
MPS  to  generate  a  picture  of  very  short  patterns  and  changes  in  atmospheric 
variables  in  the  lower  troposphere  can  lead  to  a  better  understanding  of  the 
atmosphere  and  better  modeling  at  smaller  scales. 

Future  plans  for  the  MPS  include  development  of  a  much  more  compact  ver¬ 
sion  with  improved  mobility  and  capability. 

The  MPS  has  operated  successfully  in  different  climates,  including  operations 
during  the  Los  Angeles  Free  Radical  Experiment  (LAFRE)  in  Claremont,  CA, 
and  tests  at  White  Sands  Missile  Range  (WSMR),  NM,  Erie,  CO,  Ft.  Sill,  OK, 
and  Wallops  Island,  VA. 
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1.  Introduction 


1.1  Background 

Ground-based  systems  currently  in  use  for  operational  measurement  of  atmos¬ 
pheric  profiles  rely  heavily  on  balloon-borne  rawinsondes.  The  time  between  bal¬ 
loon  launches  may  be  as  little  as  1  to  4  h  during  field  experiments,  but  in  normal 
operations  launches  occur  every  12  hours.  The  National  Oceanic  and 
Atmospheric  Administration’s  (NOAA)  National  Profiler  Network  (NPN),  con¬ 
sisting  of  404  MHz  radars  at  fixed  sites  mostly  in  the  central  United  States,  can 
provide  wind  profiles  every  hour.  [1,  2,  3,  4]  Some  sites  are  also  equipped  with 
Radio  Acoustic  Sounding  Systems  (RASS),  providing  profiles  of  virtual  temper¬ 
ature  (Ty).  Meteorological  satellite  sounders  using  either  infi-ared  or  microwave 
wavelengths  provide  a  means  of  obtaining  atmospheric  soundings  on  a  routine 
basis  for  regions  where  surface  and  upper-air  stations  are  absent.  However,  for 
mesoscale  areas  over  land,  satellite  sounder  data  may  have  horizontal  and  vertical 
resolutions  that  are  too  coarse  for  certain  applications,  especially  for  the  lower  tro¬ 
posphere.  Orlanski  defines  various  mesoscale  size  ranges.  [5] 

In  the  lowest  1  or  2  km  of  the  troposphere  over  land,  satellite  temperature  (T) 
soundings  without  ancillary  data  generally  have  errors  of  as  much  as  5  to  8  °K. 
[6,  7,  8]  For  derived  variables  such  as  wind  velocity  the  margin  for  error  may 
be  even  greater.  [9]  Vertical  and  horizontal  resolutions  are  typically  about  3  to 
5  km  and  30  to  200  km,  respectively,  depending  on  whether  infi-ared  or 
microwave  sounders  are  used  and  on  the  amount  of  spatial  averaging  for  noise 
reduction.  [10,  11,  12] 

1.2  The  Mobile  Profiler  System 

The  Mobile  Profiling  System  (MPS),  currently  being  developed  by  the  U.  S. 
Army  Research  Laboratory  (ARL)  and  NOAA’s  Environmental  Technology 
Laboratory  (ETL),  can  provide  soundings  in  the  troposphere  as  often  as  once 
every  3  min.  [13]  The  radar  wind  profiler  operating  at  924  MHz  can  provide 
wind  profiles  with  a  vertical  resolution  of  100  m  up  to  an  average  height  of  3  to 
5  km,  depending  on  atmospheric  conditions.  Under  certain  atmospheric  condi¬ 
tions  (i.e.,  moist  and  turbulent)  heights  over  6  km  are  possible.  The  RASS  can 
produce  soundings  of  Ty  up  to  0.8  to  1 .6  km  at  a  vertical  resolution  of  about 
100  m  (again  depending  on  atmospheric  conditions).  A  microwave  radiometer 
operating  in  the  oxygen  band  fi-om  50  to  60  GHz  is  able  to  produce  useful  Ty 
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profiles  to  altitudes  of  3  to  5  km.  A  second  radiometer  produces  estimates  of  total 
water  content  (vapor  and  liquid).  A  new  radiometer  currently  under  evaluation 
will  replace  both  older  radiometers  in  a  package  smaller  than  either. 

The  MPS  receives  direct  read-out  data  from  the  NOAA’s  series  of  polar-orbiting 
satellites.  The  satellite  receiver  and  processor  system  is  being  upgraded  to  a 
smaller  but  more  capable  version  that  will  be  able  to  obtain  direct  read-out  data 
from  Defense  Meteorological  Satellite  Program  (DMSP)  satellites. 

The  MPS  has  certain  elements  in  common  with  fixed-site  systems  described  by 
Parsons  et  al.  and  Stokes  and  Schwartz,  but  it  also  has  a  number  of  additional 
features.  [14,  15]  These  additions  include  software  for  processing  and  quality¬ 
controlling  data  from  the  ground-based  sensors,  and  for  combining  satellite 
soundings  with  ground-based  profiles  in  near  real  time.  Wolfe  et  al.  provide 
details  about  the  MPS  as  it  was  configured  and  operated  during  the  Los  Angeles 
Free  Radical  Experiment  (LAFRE)  in  Claremont,  CA,  and  present  examples  of 
the  various  data  processing  and  output  available.  [13]  Cogan  and  Izaguirre  pre¬ 
sent  additional  samples  of  output  and  give  preliminary  quantitative  results.  [16] 
Figure  1,  from  Wolfe  et  al.,  shows  the  primary  sensors  in  the  MPS  as  configured 
during  the  LAFRE.  [13]  Table  1  presents  certain  MPS  instmment  characteris¬ 
tics  and  compares  the  LAFRE  configuration  with  the  MPS  with  recent  and  ongo¬ 
ing  changes. 


MOBILE  PROFILER  SYSTEM 


11 


Figure  1.  Layout  of  primary  components  of  the  MPS,  from  Wolfe  et  al.  (1995). 


Table  1.  MPS  sensors  and  characteristics.  Initial  configuration  com¬ 
pared  with  that  with  recent  and  ongoing  upgrades.  New  ver¬ 
sion  has  smaller,  more  robust  shelter  for  processors  and  some 
instruments 


System 

Initial  MPS  as  during  LAFRE 

MPS  with  new/ongoing  upgrades 

Radar  wind 
profiler 

RASS 

924  MHz,  phased  array 

-2000  Hz,  4  external  sources 

Combined  phased  array  wind  radar 
(924  MHz)  and  RASS  (-2000  Hz), 
120  transducers  in  RASS 

MW 

Radiometer 

T:  50-60  GHz  (O2);  PW,  LW: 

22.7,  3 1.4  GHz,  2  radiometers 
(T,  PW/LW) 

One  radiometer  opertating  in  same 
frequency  bands  for  some  vari¬ 
ables,  smaller  than  either  of  older 
radiometers 

Satellite  receiver 
and  processor 

HRPT  from  NOAA  satellites 
(soundings  and  imagery) 

Upgraded,  smaller  system:  HRPT 
from  NOAA  and  direct  read-out 
DMSP 

Portable  surface 
station 

Standard  meteorological  variables, 
on  trailer-mounted  mast 

No  change 

Weather  map 
receiver 

Receive  maps  and  other  data  via 
GOES 

Eliminated  GOES  link,  maps  and 
other  weather  data  via  Internet 

GPS  receiver 

Provides  site  location 

No  change 

1.3  Overview 


The  merging  method  for  combining  ground-based  and  satellite  profiles 
described  in  this  report  is  a  revised  version  of  the  technique  described  in  Cogan 
and  Izaguirre.  [17]  This  method  may  be  used  for  T^,  pressure,  wind  velocity, 

and  other  meteorological  variables.  Even  though  the  microwave  radiometer 
component  of  the  MPS  uses  a  statistical  method  for  retrieval  of  temperature  pro¬ 
files  and  moisture  parameters  that  requires  a  priori  data,  the  merging  algorithm 
itself  may  be  applied  wherever  the  MPS  is  located  (i.e.,  it  is  not  site-specific  and 
requires  no  a  priori  information).  However,  the  user  may  alter  certain  software 
parameters  (e.g.,  output  layer  thickness  or  maximum  distance  fi’om  the  MPS  site 
for  acceptance  of  satellite  profiles).  Current  statistical  techniques  for  merging 
ground-based  and  satellite  profiles  of  T  (or  Ty)  reported  in  the  literature  are  site- 
specific  in  that  statistical  coefficients  are  computed  using  a  large  set  of  a  priori 
data  normally  gathered  for  a  long  series  of  rawinsonde  soundings  from  a  partic¬ 
ular  location.  [18, 19, 20]  In  this  report  we  discuss  characteristics  of  instruments 
of  the  type  employed  in  the  MPS  and  describe  the  method  for  merging  satellite 
and  ground-based  profiles  into  a  combined  sounding.  We  present  results  of  field 
tests  of  both  combined  soundings  and  profiles  from  individual  ground-based, 
sensors.  Comparisons  with  rawinsondes  and  radar-tracked  pibals  give  an  idea  of 
how  these  systems  compare  with  more  traditional  sounding  systems. 
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2.  Sensor  Characteristics 


Before  evaluating  the  merging  algorithm  and  accuracies  of  the  component  sen¬ 
sors,  the  accuracies  and  other  relevant  measurement  parameters  of  similar  data 
sources  from  the  formal  literature  should  be  examined.  Of  particular  interest 
are  the  satellite  sounders  and  ground-based  radar  profilers,  RASS,  and 
microwave  radiometers  of  the  type  used  in  the  MPS.  Table  2  summarizes 
selected  information  extracted  from  Moran  and  Strauch,  Flowers  et  al., 
Okrasinski  and  Olsen,  Weber  and  Wuertz,  Weber  et  al..  Reale,  May  et  al.,  Le 
Marshall,  Franklin  and  Lord,  and  Jedlovec.  [21,  22,  23,  24,  3,  8,  25,  7,  9,  6] 
These  values  may  be  compared  with  similar  data  for  common  rawinsonde  sys¬ 
tems  (Fisher  et  al.),  also  presented  in  table  2.  [26] 


Table  2.  Capabilities  of  several  remote  sensing  and  rawindsonde  systems. 
The  values  shown  represent  averages  for  most  current  systems  in  each  group. 


System 


Radar 
Profiler 
915/924  MHz 


RASS 

Microwave 

Radiometer 


Variable 


wind  speed 

wind 

direction 

v&tual 
temperature 
temp,  or 
virtual  temp. 


Accuracy 

±  1.5to3ms“* 

±  10  to  15° 

±lto2°K 

±ltD3°K 


Capabilities 

Temporal 

resolution 


3  to  6  min 
(Consensus 
methods  may 
need  up  to 
30  min) 


3  min 

“h(2  Sat)“ 
4  h  (3  Sat) 
(1  to  2  h  near 
poles) 

1  to  3  h 
depending  on 
maximum 
height 


^  ^  -  temperature  ±  2  to  2.5°  K 

Satelhte  Sounder 

(TOYS,  SSM/T-1)  wind  speed  ±  4  to  14  ms-i 
wind  direction  ±  10  to  30° 

temperature  ±  0,5  to  1°  K 
Rawindsonde  wind  speed  ±  0.5  to  2  ms-* 
wind  direction  ±  5  to  10° 


Vertical  Vertical 

resolution  range 


I  100  to  5000  m 

100  m 


100-1600  m 


variable  up  to  10  km 


3  to  5  km  2  to  40  km 


point  vahie 

60to600m  Surfeceto30km 
layers 


From  table  2  we  see  that  ground-based  radar  profilers  and  RASS  have  accura¬ 
cies  that  approach  those  of  some  operational  rawinsonde  systems.  Flowers  et 
al.  note  that  the  RASS  performance  may  be  reduced  under  certain  conditions 
(e.g.,  strong  near-surface  winds),  and  the  extreme  cases  they  reported  are  not 
included  in  the  table.  [22]  All  the  systems  in  table  2  produce  measurements  of 
a  volume  dr  layer  average,  except  rawinsonde  values  of  T.  A  radar  profiler 
obtains  a  mean  wind  over  a  volume,  with  a  horizontal  scale  on  the  order  of  tens 


15 


to  thousands  of  meters  depending  on  altitude  from  the  radar,  beam  elevation 
angle,  and  beamwidth.  For  a  system  such  as  404,  449,  or  50  MHz  radars,  the 
volume  represented  by  the  three  (or  five)  beams  used  to  measure  the  horizontal 
and  vertical  wind  may  have  a  horizontal  diameter  exceeding  10  km  near  15  or 
20  km  altitude.  The  radar  processing  algorithms  developed  by  ETL  for  the 
MPS  (Wolfe  et  al.)  represent  an  improvement  over  the  standard  consensus  tech¬ 
niques.  [13]  These  algorithms  allow  for  higher  data  quality  at  faster  data  rates 
than  previously  possible.  Errors  in  wind  measurements  caused  by  migratory 
birds  (Wilczak  et  al.)  have  been  addressed  through  a  new  spectral  averaging 
method  (Merritt).  [27,  28]  Although  limited  when  bird  densities  are  high,  this 
algorithm  attempts  to  identify  interference  in  the  spectral  data  prior  to  averag¬ 
ing,  thereby  retaining  more  useful  wind  information. 


Vertical  profiles  of  Ty  can  be  inferred  from  measurements  of  microwave  bright¬ 
ness  temperatures.  For  surface-based  measurements,  in  the  20  to  60  GHz 
region,  the  measured  brightness  temperatures  approximately  satisfy  the  follow¬ 
ing  equation: 


T(z)a^(z)exp[  -j  \(z  )dz  ''[dz + T^yCxp  -j" 


a^(z)dz 


where  T^y  is  the  downwelling  microwave  brightness  temperature  at  frequency 

v;  Tfz;  is  temperature  at  height  z;  a/z)  is  the  absorption  coefficient;  T”^y  is  the 
downwelling  cosmic  microwave  background  brightness  temperature  above  the 
atmosphere. 

Inferring  atmospheric  temperature  structure  from  microwave  brightness  tem¬ 
perature  measurements  thus  becomes  the  problem  of  solving  (inverting)  equa¬ 
tion  (1)  to  find  T(zJ.  A  database  of  past  radiosonde  observations  has  been  used 
to  calculate  corresponding  received  radiances  at  our  operating  frequencies.  Our 
temperature  profiles  are  calculated  with  regression  coefficients  computed  from 
the  database  radiances  and  corresponding  radiosonde  observations.  The  absorp¬ 
tion  coefficient  cCyfzJ  in  the  frequency  region  of  interest  is  due  mainly  to  oxy¬ 
gen  lines  (50  to  60  GHz),  the  water  vapor  line  at  22.235  GHz,  and  a  liquid  water 
continuum  measured  near  3 1 .4  GHz.  Oxygen  is  well  mixed  in  the  atmosphere, 
so  a  good  a  priori  estimate  of  its  contribution  is  possible.  In  our  scheme,  water 
vapor  and  liquid  water  are  independently  measured  by  radiometric  channels 
near  22.235  and  31.4  GHz.  Vertical  resolution  is  about  30  m  at  heights  below 
1  km,  increasing  to  more  than  1  km  around  10  km. 


Satellite  sounders  measure  radiances  that  are  converted  to  temperatures  that 
represent  means  for  large  volumes  of  atmosphere,  according  to  the  frequency- 
dependent  weighting  functions  and  horizontal  field  of  view.  Temporal  resolu¬ 
tions  in  table  2  are  average  values.  Generally,  satellite  values  are  valid  for  alti¬ 
tudes  >  2  or  3  km  above  a  land  surface  and  not  near  the  tropopause  (e.g.,  T  to 
+  5  to  8  °K  near  the  surface,  and  +  3  to  5  °K  near  the  tropopause).  A  tempera¬ 
ture  may  represent  the  mean  over  a  vertical  extent  of  3  to  5  km  and  over  a  hor¬ 
izontal  area  of  tens  to  over  100  km  diameter  (assuming  a  circular  area).  Wind 
velocity  is  derived  from  the  satellite  temperature  profile  using  geostrophic,  gra¬ 
dient,  or  thermal  wind  equations  (Franklin  and  Lord).  [9]  A  rawinsonde 
acquires  a  mean-layer  wind  velocity  along  its  path,  over  a  period  of  perhaps 
1  min  or  less  (most  rawinsondes  rise  at  about  5  ms‘^).  Its  temperature  mea¬ 
surements  may  be  considered  as  point  values.  Temporal  resolution  of  <  1  h  may 
be  achieved  for  rawinsondes  if  there  is  a  capability  for  multiple  transmitter  fre¬ 
quencies.  Table  2  provides  an  indicator  of  the  relative  “quality”  of  data  from 
the  listed  sensors  inclusive  of  inherent  differences  between  them  due  to  mea¬ 
surement  and  processing  methods.  However,  satellite  sounders  are  the  only 
means  of  obtaining  large-area  or  global  coverage,  while  radar  profilers  with 
RASS  and  microwave  radiometers  have  the  best  spatial  and  temporal  resolution. 
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3.  Merging  Algorithms 


The  algorithms  described  in  this  section  combine  profiles  firom  a  suite  of 
ground-based  systems  and  satellite  sounders.  Radar  profiler,  RASS,  and 
microwave  radiometer  provide  data  below  the  lowest  satellite  level.  Currently, 
RASS  values  of  Ty  are  used  up  to  the  highest  RASS  data  level.  Above  that 

height  radiometer  values  are  used,  when  available,  up  to  the  maximum  height 
of  “useful”  data  (nominally  about  3.5  km  based  on  early  test  results).  The  com¬ 
bined  RASS  and  radiometer  profile  is  the  ground-based  profile  of  Ty.  Where 
ground-based  and  satellite  profiles  overlap,  the  satellite  data  are  weighted  in 
accordance  with  the  spatial  and  temporal  separation  of  the  sounding  from  the 
ground-based  profiles.  The  radar  wind  profiler,  RASS,  and  radiometer  combi¬ 
nation  will  hereafter  be  referred  to  collectively  as  the  profiler.  The  spatial 
weighting  function  has  an  elliptical  form: 


w=\-x'^ia'^-y^ib^ 


(2) 


where  W  is  weight,  A  and  B  are  the  semi-major  and  semi-minor  axes  of  the 
assumed  elliptical  area  represented  by  the  satellite  sounding  (>  the  horizontal 
resolution  of  a  single  satellite  profile),  and  X  and  Y  are  the  distances  along  those 
axes  fi’om  the  profiler  to  the  center  of  the  satellite  sounding  footprint.  For 
polar-orbiting  satellites  a  circular  area  may  be  sufficient  away  fi-om  the  edges 
of  the  swath,  say  within  500  or  600  km  of  the  subsatellite  track.  Some  misreg¬ 
istration  with  height  can  occur  as  nadir  angle  increases,  especially  towards  the 
edges  of  the  swath.  The  temporal  weighting  function  has  an  initial  period  when 
the  two  sources  of  data  have  equal  weight  (e.g.,  15  min),  followed  by  a  period 
of  linear  decrease  to  some  time  (usually  3  to  6  h)  when  the  satellite  data  are 
ignored  (temporal  weight  =  0).  The  final  weight  given  to  the  satellite  sounding 
is  the  product  of  the  spatial  and  temporal  weights  times  an  accuracy  ratio  (R). 
This  ratio  relates  typical  accuracies  (table  2)  of  the  radar  wind  profiler,  RASS 
(Ty),  and  radiometer  (T  or  Ty  estimate)  to  those  of  the  satellite  sounder.  For 

current  instruments,  R  decreases  the  weight  given  to  the  satellite  data.  The  user 
may  alter  these  parameters. 
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The  satellite  and  ground-based  profiles  may  overlap,  or  a  gap  may  occur 
between  them.  When  the  profiles  overlap,  the  satellite  data  are  interpolated  to 
profiler  heights  in  the  overlap  region.  The  equation  for  combining  the  two  sets 
of  data  for  those  heights  has  the  following  form: 

Q=ff^(Qp^Q,y2^(i-fV)Q^  (3) 


where  0  is  a  variable  at  some  height,  fV  is  here  the  combined  temporal  and  spa¬ 
tial  weight  times  the  accuracy  ratio  R,  and  subscripts  ^  and  ^  refer  to  profiler  and 
satellite,  respectively. 

Where  the  data  do  not  overlap,  a  gap  exists  between  the  highest  altitude  of  a 
ground-based  profile  and  the  lowest  height  of  the  satellite  data.  In  this  case  the 
satellite  data  are  extrapolated  down  to  the  maximum  altitude  of  the  ground- 
based  profile.  Above  the  highest  altitude  of  the  profiler  data  for  either  T^  or 
wind  velocity,  the  satellite  value  is  adjusted  according  to  a  scheme  described  in 
Cogan  and  Izaguirre.  [17]  It  is  based  on  the  difference  between  satellite  (actu¬ 
al  or  extrapolated)  and  profiler  values  at  the  maximum  profiler  height  for  the 
particular  variable.  The  adjustment  or  correction  is  reduced  in  magnitude 
through  multiplication  by  an  adjustment  parameter  (a)  at  successive  heights  up 
to  a  preset  number  of  satellite  levels  (i.e.,  D,  =  afDj.j],  where  D  is  the  differ¬ 
ence  value,  and  ,■  is  a  satellite  sounding  level).  Normally  data  are  adjusted  for 

three  to  five  satellite  levels  above  the  highest  profiler  level.  The  adjustment 

parameter  a  and  the  number  of  levels  may  be  altered  by  the  user.  If  the  gap 
between  the  highest  profiler  level  and  the  satellite  level  immediately  above 
exceeds  a  preset  value  (e.g.,  2  or  3  km)  the  algorithm  skips  the  extrapolation 
routine  and  does  not  adjust  the  satellite  data.  Each  satellite  T  profile  is  con¬ 
verted  to  Ty  using  retrieved  dew  points,  if  available;  otherwise  the  program  uses 

a  rough  estimate  based  on  the  surface  value  of  humidity.  Alternatively,  a  profile 
based  on,  say,  regional  climatology  could  be  used.  At  the  heights  of  satellite 
data  used  here,  >  2.5  km  above  ground  level  (AGL),  Ty  often  is  within  1  °K  of  T. 
No  conversion  takes  place  if  T  or  dewpoint  is  less  than  233  °K,  or  z  >  10  km. 

Figure  2  shows  plots  of  combined  soundings  of  wind  velocity  for  the  period  170 
to  2000  UTC,  July  29,  1994,  at  White  Sands  Missile  Range  (WSMR),  NM. 
Each  wind  profile  is  a  1 5  min  average  ending  at  the  time  when  the  profile  is  plot¬ 
ted.  Wind  barbs  plotted  near  4.8,  6.4,  and  8.3  km  are  satellite-derived  thermal 
winds  modified  according  to  the  aforementioned  method. 
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4.  Results  and  Comparisons 


4.1  Los  Angeles  Free  Radical  Experiment 

Personnel  from  ARL  and  ETL  participated  in  the  Los  Angeles  Free  Radical 
Experiment  (LAFRE),  using  the  MPS  to  obtain  detailed  sounding  data  for  the 
primary  sponsor,  the  California  Air  Resources  Board.  These  data  also  served 
to  test  the  system  and  algorithms.  The  MPS  operated  almost  continuously  from 
August  28  through  September  23,  1993.  During  this  period  a  MARWIN  and 
Cross-chain  Loran  Atmospheric  Sounding  System  (CLASS)  were  operated 
simultaneously  from  the  MPS  trailer.  The  ability  to  operate  rawinsondes  from 
the  MPS  allowed  for  detailed  intercomparisons.  Near  real-time  graphical  and 
statistical  comparisons  were  possible  using  software  developed  by  ETL  and 
ARL  (Wolfe  et  al.;  Cogan  and  Izagirrre).  [13,  16]  The  microwave  radiometer 
(temperature  only)  operated  only  during  part  of  the  LAFRE,  and  the  merged 
profiles  discussed  here  do  not  include  microwave  data. 

From  August  28  through  September  11,  1993,  the  Los  Angeles  basin  was  under 
a  strong  upper  ridge,  and  at  times  a  closed  high  pressure  area  from  the  surface 
through  300  hPa.  The  marine  boundary  layer  was  consistently  capped  by  one  or 
more  inversions.  Wolfe  et  al.,  and  Cogan  and  Izaguirre  present  charts  that  show 
wind  velocities  from  the  radar  profiler  for  typical  days  during  this  early  part  of 
the  LAFRE,  depicting  light  and  often  variable  winds.  [13,  16]  Combining  these 
profiles  with  the  nearest  good  satellite  sounding,  sometimes  as  much  as  300  km 
distant,  led  to  a  “worst  case”  situation  on  several  days.  Atmospheric  conditions, 
especially  wind  velocity,  are  often  quite  different  200  or  300  km  to  either  side  of 
a  strong  ridge. 

Figures  3  through  6  show  merged  profiles  from  the  LAFRE  data  compared  with 
rawinsonde  data  for  0.1  km  averaged  layers.  Figure  3  contains  up  to  36  poten¬ 
tial  comparisons  during  September  7  through  11,  compared  with  a  maximum  of 
12  during  September  17  and  20  through  23  for  figure  4.  During  the  second  peri¬ 
od,  comparisons  were  obtained  on  September  1  and  September  20  through  23. 
Figures  5  and  6  show  the  RASS  values  of  figures  3  and  4,  respectively,  in  an 
expanded  scale  for  heights  <  2  km.  For  these  comparisons  RASS  T^  values  were 

corrected  for  vertical  velocity  in  a  manner  similar  to  that  employed  by  Moran 
and  Strauch.  [21]  These  results  suggest  that  the  greatest  differences  may  occur 
within  any  layer  of  atmosphere,  and  are  not  uniformly  distributed  with  height. 
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Figure  3.  Mean  (solid  line)  and  standard  deviation  (dashed  line)  of  (°K) 


differences  (MPS  -  rawindsonde)  for  September  7  through  11, 
1993.  Heights  are  AGL. 
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Figure  6.  Same  as  figure  4,  but  for  RASS  only. 
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Figure  3  shows  a  fairly  typical  pattern  (e.g,,  accuracy  degraded  near  the 

tropopause)  of  standard  deviation  (Oj^)  of  differences  at  satellite  levels  with 
values  somewhat  lower  than  previously  published  (Le  Marshall;  Reale).  [7,  8] 
The  merging  of  satellite  and  RASS  profiles  greatly  reduced  the  differences  in 
the  lowest  few  kilometers  relative  to  satellite  soundings  alone  (about  0.5  to  2  °K 
versus  5  to  8  °K).  Figures  3  and  5  show  that  the  magnitude  of  mean  differences 
and  Oy- values  for  the  RASS  are  fairly  small  (<  0.5  and  <  0.9  °K,  respectively, 
for  0.1  <  z  <  1.3  km),  except  at  the  highest  RASS  levels  (>  1 .3  km)  and  near  the 
surface  (<  0.1  km). 

Figure  4  presents  a  less  common  distribution  of  values  with  height  (z)  at  satel¬ 
lite  levels.  The  O'p  of  Ty  differences  decreases  near  the  tropopause,  and  the 
mean  of  the  differences  is  high  at  lower  satellite  levels  in  spite  of  a  reduction  in 
magnitude  as  a  result  of  the  merging  process.  In  figures  4  and  6  the  RASS  val¬ 
ues  of  Of  are  higher  at  z  <  0.5  km  and  at  z  >  1.3  km.  Magnitude  of  the  mean 
differences  for  RASS  is  high  at  some  lower  layers,  and  at  1 .6  km.  These  results 
are  not  unexpected.  In  the  lower  layers,  surface  heating  and  cooling  cause 
greater  variability.  The  higher  layers,  near  the  maximum  range  for  RASS,  often 
coincided  with  the  height  of  the  marine  inversion.  Variability  in  height  and 
strength  of  this  inversion,  and  its  relationship  to  location  and  thickness  of  com¬ 
bined  sounding  layers,  are  significant  factors  when  comparing  with  rawinsonde 
soundings. 

The  lower  Oj  near  the  tropopause  during  the  second  period  may  have  arisen  as 
a  result  of  a  weaker  lapse  rate,  or  little  change  in  it,  and  possibly  also  as  a 
result  of  a  lower  horizontal  temperature  gradient  near  the  tropopause  over 
the  area  in  and  around  the  Los  Angeles  basin  during  most  or  all  of  the  sec¬ 
ond  period.  The  former  would  tend  to  reduce  inaccuracies  in  the  satellite  data 
due  to  the  inherent  smoothing  of  vertical  temperature  gradients  and  errors 
induced  by  incorrect  height  attributed  to  a  satellite  value.  The  latter  would  tend 
to  lower  temperature  differences  between  satellite  and  rawinsonde  that  occur 
because  of  horizontal  distance  between  satellite  sounder  field  of  view  and 
rawinsondJe  location.  An  investigation  of  sounding  data  fi’om  three  stations 
near  the  Los  Angeles  basin  indicated  that  the  magnitude  of  the  lapse  rate  fi’om 
about  11  to  14  km  (below  the  tropopause)  was  about  1.2  °Kkm’i  to  about  1.7 
°Kkm‘i  smaller  during  the  comparison  days  of  the  second  period  (average) 
depending  on  the  station.  Also,  the  tropopause  was  about  2  to  3  °K  warmer  dur¬ 
ing  the  second  period. 

Upper-air  maps  for  200  hPa  (12.1  km  <  z  <  12.5  km),  near  the  top  of  the  com¬ 
bined  soundings  from  the  LAFRE,  showed  a  more  zonal  distribution  of  tern- 


perature  during  the  second  period  than  during  the  first.  The  mean  and  standard 
deviation  of  the  magnitude  of  the  horizontal  lapse  rate  within  300  km  of  the  test 
site  showed  almost  no  change  for  the  east-west  gradient.  Small  increases  in 
magnitude  of  mean  and  standard  deviation  of  about  0.1  to  0.2  °K/100  km  each 
appeared  for  the  north-south  gradient  from  the  first  to  second  period.  An 
increase  in  the  north-south  gradient  would  be  expected  as  the  distribution  of 
temperature  became  more  zonal.  However,  the  smoothing  of  horizontal  tem¬ 
perature  gradients  in  satellite  data  (Jedlovec)  would  tend  to  further  reduce  the 
effect,  if  noticeable,  of  those  small  changes.  [6]  While  the  effect  on  tempera¬ 
ture  differences  between  satellite  and  rawinsonde  would  have  been  minimal 
from  the  changes  in  the  horizontal  gradient  noted  above  (very  small  increase  in 
Oj),  the  more  significant  decrease  in  vertical  temperature  gradient  and  the 
warmer  tropopause  during  the  second  period  (relatively  large  decrease  in  Oj) 
may  at  least  partly  explain  the  reduction  in  Oj  in  the  few  kilometers  below  the 
tropopause. 

Table  3  shows  means  and  standard  deviations  of  wind  speed  differences  in  meters 
per  second  for  the  radar  profiler  and  satellite  (adjusted  at  the  lowest  three  satellite 
data  levels)  relative  to  rawinsonde.  The  maximum  number  of  data  comparisons  by 
layer  are  the  same  as  for  the  T^  comparisons  presented  above. 


Table  3.  Means  and  standard  deviations  of  wind  speed  differences  (ms'i) 
for  0.3  km  layers  (indicated  sensor  -  rawinsonde)  obtained 
during  the  LAFREE. 


September  7-11, 1993 

September  17  and  20-23,  1993 

Radar  profiler 

Satellite 

Radar  profiler 

Satellite 

Mean 

0.75 

10.84 

1.53 

8.61 

Std.  deviation 

1.88 

8.60 

2.75 

2.88 

No.  of  layers 

14 

6 

11 

6 

At  times  the  rawinsonde  data  may  contain  serious  errors.  Fisher  et  al.  present 
information  on  the  average  errors  found  in  several  types  of  rawinsonde  systems. 
[26]  To  gain  an  idea  of  the  quality  of  the  rawinsonde  data  from  the  LAFRE, 
soundings  were  compared  from  two  similar  systems  (MARWIN  and  CLASS) 
receiving  data  from  one  sonde.  Differences  in  T^  from  comparisons  using  a  sin¬ 
gle  sonde  averaged  approximately  ±  0.2  to  0.4  °K,  with  maximum  differences 
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of  about  ±  1  °K.  Figure  7  compares  wind  speed  differences  between  the  two 
systems.  The  periodic  pattern  is  consistent  with  other  data  examined  to  date. 
The  large  differences  near  and  above  3  km  are  on  the  high  side,  but  values 
around  ±  1  ms'*  are  not  uncommon.  Cogan  and  Izaguirre  presented  data  show¬ 
ing  a  few  wind  direction  variations  of  >  90°  in  one  case  during  the  LAFRE, 
although  wind  direction  differences  for  most  heights  in  data  examined  to  date 
were  <  10°.  [16]  This  type  of  comparison  suggests  that  differences  in  profiler 
wind  speed  and  direction  of  around  ±  1  ms*’  and  10°,  respectively,  relative  to 
rawinsonde  may  be  close  to  the  “best”  one  could  expect.  A  possible  partial 
explanation  for  the  wind  speed  differences  is  that  the  MARWTN  software  has 
more  extensive  built-in  checks  and  smooths  the  data  somewhat.  Nevertheless, 
caution  must  be  taken  when  using  a  rawinsonde  sounding  as  a  standard,  espe¬ 
cially  in  light  winds.  The  user  should  make  sure  each  sounding  contains  valid 
data,  and  apply  appropriate  quality  controls. 

4.2  Wallops  Island 

Field  tests  at  the  NASA  Wallops  Flight  Facility  (WFF)  on  Wallops  Island,  VA, 
provided  the  opportunity  to  compare  MPS  wind  profiles  for  the  lowest  1 .9  km 
with  wind  profiles  obtained  jfrom  radar-tracked  pibal  balloons.  An  unusual 
aspect  of  this  experiment  as  compared  with  other  comparison  studies  such  as 
the  LAFRE  and  earlier  work  (e.g.,  Weber  and  Wuertz)  was  the  ability  to  exam¬ 
ine  the  background  wind  variability  at  the  same  time  as  the  comparisons.  [24] 
During  the  week  of  July  17  through  21,  1995,  for  morning  and  afternoon  peri¬ 
ods  lasting  about  1  to  1  Vi  h,  two  pibals  were  launched  about  3  min  apart  every 
1 5  min  (four  to  five  “pairs”  each  period).  The  MPS  operated  continuously  dur¬ 
ing  these  periods  producing  wind  profiles  every  3  min.  For  each  pibal  pair, 
comparisons  were  made  between  the  3  min  MPS  profile  just  prior  to  the  second 
pibal  and  the  second  pibal,  and  between  the  first  and  second  pibal.  Surface  val¬ 
ues  shown  were  taken  from  the  WFF  and  MPS  surface  sensors.  The  site  of  the 
experiment  was  about  0.2  km  west  of  the  ocean,  with  the  MPS  located  <  50  m 
east  from  the  pibal  launch  site.  Figures  8  through  11  show  the  means  and  stan¬ 
dard  deviations  of  the  wind  speed  and  direction  differences  between  MPS  and 
second  pibal,  and  pibal  pairs  for  100  m  layers,  during  July  20, 1995  (nine  pairs), 
and  July  21,1995  (five  pairs). 

The  MPS  versus  pibal  comparison  for  July  1 8  (not  shown)  showed  significantly 
greater  differences  in  wind  direction  than  for  the  other  two  comparison  days 
(July  20  and  21).  The  pibal  versus  pibal  comparison  for  that  day  also  showed 
somewhat  larger  differences  in  wind  direction  relative  to  those  for  the  other  days. 


Figure  7.  Mean  wind  speed  differences  (ms**)  between  CLASS  and 

MARWIN  systems  from  LAFRE  data.  Data  gathered  by  one 
sonde.  Heights  are  AGL. 


Figure  8.  Mean  wind  speed  differences  (ms’*)  from  WFF  data.  Heights 

are  AGL.  Bold  curves  represent  MPS  versus  pibal,  lighter  ones 
represent  pibal  versus  pibal.  Left  graph  (square  data  points)  and 
right  graph  (circles)  are  plots  for  July  20  and  21,  respectively. 
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Figure  9.  Same  as  flgure  8,  except  plotted  for  standard  deviation  of  wind 
speed  differences  (ms**). 
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Figure  11.  Same  as  figure  8,  except  plotted  for  standard  deviation  of 
wind  direction  differences  (°). 
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The  standard  deviation  of  wind  direction  differences  exceeded  10°  for  pibal 
versus  pibal  for  heights  (z)  <  0.3  km  and  z  =  1.1  km,  and  for  the  MPS  versus 
pibal  at  z  <  0.8  km  and  z  =  1 .2  km.  The  magnitude  of  mean  differences  in  wind 
direction  between  MPS  and  pibal  was  >  20°  at  z  =  0.6  km,  0.8  <  z  ^  1.2  km, 
and  1.5  <  z  <  1.8  km  (maximum  of  about  25°  at  1.6  km).  Magnitude  of  mean 
and  standard  deviation  of  wind  speed  differences  was  somewhat  larger  at  z  < 
0.6  km  for  the  MPS  versus  pibal.  Also,  on  this  day  the  pibals  traveled  eastward, 
passing  over  the  ocean  within  1  min  after  launch.  These  larger  differences 
were,  therefore,  not  unexpected  since  the  pibals  drifted  over  the  ocean  after 
reaching  200  or  300  m  in  altitude,  leaving  the  highly  convective  conditions  that 
existed  over  the  land.  Later  in  the  afternoon  small  but  intense  thunderstorms 
passed  through  fi-om  the  west,  forcing  the  test  to  be  canceled  before  1500  EDT 
(1900  UTC)  due  to  the  danger  of  lightning  strikes. 

The  largest  direction  difference  between  the  MPS  and  pibal  was  at  the  “surface” 
(about  5  m  AGL)  for  the  latter  two  days  and  at  0.1  km  on  July  18.  Both  sys¬ 
tems  relied  on  surface  stations  separated  by  about  10  m  horizontally  and  1  to 
2  m  vertically  (the  WFF  anemometer  was  higher).  The  location,  only  about 
200  m  fi’om  the  ocean,  and  the  mix  of  land  and  water  surfaces  near  the  launch 
site  may  account  for  much  of  the  observed  direction  differences  in  the  lowest 
0.1  to  0.2  km.  The  balloons  drifted  off  roughly  to  the  northwest  except  on 
July  1 8  when,  soon  after  turning  toward  the  east  to  northeast,  they  passed  over 
the  northern  half  of  the  island  and  then  out  over  the  water.  Since  the  ascent  rate 
of  the  pibals  was  about  5  ms‘*  and  the  average  wind  speed  for  most  of  the  test 
periods  was  about  5  to  7  ms'^  during  much  of  each  ascent,  the  balloon  ended  up 
about  2  to  3  km  from  the  MPS  and  pibal  launch  site  by  the  time  it  reached  an 
altitude  of  2  km.  On  July  21  the  wind  speed  at  most  heights  exceeded  10  ms'^ 
causing  the  pibal  to  drift  about  4  km  by  the  time  it  rose  to  2  km. 

The  comparison  of  wind  profiling  radar  with  radar-tracked  pibals  yielded 
results  for  wind  speed  that  appear  better  than  those  shown  in  table  3.  The  aver¬ 
age  magnitude  of  mean  differences  (0.1  km  layers)  from  the  WFF  data  (fig¬ 
ure  8)  appears  similar  to  that  firom  the  LAFRE  data  for  September  7  through  11, 
1993.  However,  the  average  of  the  means  from  the  LAFRE  for  September  17 
and  20  through  23  is  considerably  larger.  The  standard  deviations  of  the  dif¬ 
ferences  for  all  layers  in  figure  9  are  less  than  the  average  standard  deviation  of 
either  period  in  table  3.  The  average  standard  deviation  between  MPS  and  pibal 
in  figure  9  is  about  0.8  ms'L 


The  accuracies  suggested  using  data  from  the  WFF  experiment  (figures  8 
through  11)  appear  better  overall  than  those  presented  in  table  2  for  radar  pro¬ 
filers,  except  for  wind  direction  on  July  1 8.  The  differences  between  values  in 
figures  8  through  11  and  tables  2  and  3,  and  variations  between  profiles  from 
pibals  launched  3  min  apart  (also  figures  8  through  11),  support  the  idea  that  at 
least  part  of  the  differences  between  radar  profiler  and  rawinsonde  wind  sound¬ 
ings  in  earlier  LAFRE  data  indeed  may  be  a  result  of  real  atmospheric  tempo¬ 
ral  and  spatial  variation.  The  data  for  1 8  July,  with  even  larger  standard  devia¬ 
tions  of  wind  direction  differences  for  both  pibal  versus  pibal  and  MPS  versus 
pibal  than  those  for  20  and  21  July,  suggests  that  these  variations  can  be  signif¬ 
icant  even  over  a  3  min  time  span. 

The  Wallops  Island  experiment  provided  an  opportunity  to  test  a  new  type  of 
surface-based  radiometer,  developed  for  ARL  by  the  OPHIR  Corporation.  This 
new  radiometer,  called  the  Next  Generation  Radiometer  (NGR),  incorporates 
several  design  advances  over  the  older  system  (the  Passive  Microwave 
Temperature  Profiler,  or  PMTP),  including  a  much  smaller  size  and  lighter 
weight,  and,  notably,  frequency  tunability  and  very  precise  software  control  of 
radiometer  frequency  in  the  oxygen  bands.  The  antenna  system  for  the  NGR  is 
based  on  an  optical  lens  that  focuses  into  a  corrugated  horn  antenna. 

The  PMTP  uses  Gunn  diodes  for  local  oscillators  for  all  measurements.  There 
are  four  such  oscillators  limiting  measurement  to  four  frequencies.  The  oscil¬ 
lators  also  suffer  from  a  tendency  to  drift  in  frequency  when  their  temperatures 
are  not  precisely  controlled,  and  require  frequent  calibration  to  check  for  drift 
due  to  mechanical  effects.  In  the  NGR,  the  oxygen-band  local  oscillator  is  a 
highly  stable  tunable  synthesizer.  This  frequency  tunability  of  NGR  makes  it 
practical  to  use  a  larger  number  of  frequencies  (for  this  experiment,  11  fre¬ 
quencies)  in  the  50-60  GHz  sensing  band. 
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In  the  Wallops  experiment,  the  performance  of  the  radiometers,  as  measured  by 
comparison  with  simultaneous  radiosonde  observations,  was  comparable  at 
lower  levels,  but  the  NGR  appeared  definitely  superior  at  the  higher  levels. 
Figure  12  shows  root  mean  square  (RMS)  differences  for  nine  comparisons 
between  simultaneous  radiosonde  observations  and  corresponding  measure¬ 
ments  with  the  current  ARL  oxygen  radiometer  system  (PMTP)  and  the  OPHIR 
radiometer  (NGR). 


Figure  12.  RMS  deviations  of  NGR  (bold  squares)  and  PMPT 

(lighter  circles)  from  concurrent  rawindsonde  observations. 


5.  Conclusions 


The  MPS  shows  promise  as  a  means  of  collecting  data  from  a  variety  of  profil¬ 
ing  instruments  and  merging  these  data  into  combined  meteorological  sound¬ 
ings  for  near  real-time  operational  applications.  The  merging  method  provides 
soundings  with  an  accuracy  in  temperature  (or  virtual  temperature)  comparable 
to  rawinsonde  soundings  or  to  other  currently  published  methods  of  combining 
ground-based  and  satellite  data.  While  the  microwave  radiometer  component 
of  the  MPS  uses  a  priori  information,  the  merging  algorithm  has  the  advantage 
of  not  being  site  specific,  and  a  priori  data  sets  are  not  required  as  in  statistical 
merging  techniques.  It  also  may  be  used  to  combine  profiles  of  other  meteoro¬ 
logical  variables.  However,  the  accuracy  of  wind  velocity  values  above  the 
maximum  radar  data  level  is  limited  by  the  errors  in  current  methods  of  deriv¬ 
ing  wind  velocity  from  satellite  sounder  data.  New  ways  to  derive  satellite 
wind  velocities  are  being  investigated.  In  the  interim,  merging  wind  data  from 
conventional  systems  (e.g.,  rawinsondes)  for  the  upper  part  of  the  sounding  may 
be  the  only  viable  alternative. 

The  combining  algorithm  is  not  limited  to  the  MPS.  It  also  can  be  applied  to 
other  suites  of  instruments  capable  of  measuring  profiles.  Existing  facilities 
where  this  algorithm  may  prove  useful  include  the  NOAA  NPN,  and  systems  at 
sites  within  or  near  airports  and  at  government  test  ranges.  The  basic  method 
may  be  applied  to  airborne  systems  as  well. 

The  data  provided  by  the  MPS  will  have  a  variety  of  civilian  and  military  appli¬ 
cations.  The  MPS  can  provide  timely  support  for  airfield  operations,  giving,  for 
example,  near  real-time  indications  of  potentially  hazardous  wind  conditions. 
As  the  Los  Angeles  Free  Radical  Experiment  showed,  this  type  of  system  can 
be  invaluable  for  pollution  studies.  The  ability  to  generate  a  picture  of  very 
short-term  flow  and  virtual  temperature  patterns  in  the  lower  troposphere  can 
lead  to  a  better  understanding  of  the  atmosphere  and  to  better  modeling  at 
smaller  scales. 
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NUC  EFFECTS  DIV 
AMSRL  SL  CO 
APG  MD  21010-5423 

US  ARMY  MATERIEL  SYST 
ANALYSIS  ACTIVITY 
AMSXY 

APG  MD  21005-5071 

ARMY  RESEARCH  LABORATORY 
AMSRL  D 

2800  POWDER  MILL  ROAD 
ADELPHl  MD  20783-1145 

ARMY  RESEARCH  LABORATORY 
AMSRL  OP  Cl  SD  TL 
2800  POWDER  MILL  ROAD 
ADELPHl  MD  20783-1 145 

ARMY  RESEARCH  LABORATORY 
AMSRL  Cl  LL 
ADELPHl  MD  20703-1197 

ARMY  RESEARCH  LABORATORY 
AMSRL  SS  SH 
ATTN  DR  SZTANKAY 
2800  POWDER  MILL  ROAD 
ADELPHl  MD  20783-1145 

ARMY  RESEARCH  LABORATORY 

AMSRL  IS 

ATTN  J  GANTT 

2800  POWDER  MILL  ROAD 

ADELPHl  MD  20783-1 197 

ARMY  RESEARCH  LABORATORY 
AMSRL 

2800  POWDER  MILL  ROAD 
ADELPHl  MD  20783-1145 


NATIONAL  SECURITY  AGCY  W21 
ATTN  DR  LONGBOTHUM 
9800  SAVAGE  ROAD 
FT  GEORGE  G  MEADE  MD  20755-6000 

US  ARMY  RSRC  OFC 
ATTN  AMXRO  GS  DR  BACH 
PO  BOX  12211 
RTP  NC  27009 

DR  JERRY  DAVIS 
NCSU 

PO  BOX  8208 
RALEIGH  NC  27650-8208 

US  ARMY  CECRL 
CECRL  GP 
ATTN  DR  DETSCH 
HANOVER  NH  03755-1290 

US  ARMY  ARDEC 
SMCAR  IMI I  BLDG  59 
DOVER  NJ  07806-5000 

ARMY  DUGWAY  PROVING  GRD 
STEDP  MT  DA  L  3 
DUGWAY  UT  84022-5000 

ARMY  DUGWAY  PROVING  GRD 
STEDP  MT  M 
ATTN  MR  BOWERS 
DUGWAY  UT  84022-5000 

DEPT  OF  THE  AIR  FORCE 

OL  A  2D  WEATHER  SQUAD  MAC 

HOLLOMAN  AFB  NM  88330-5000 

PL  WE 

KIRTLAND  AFB  NM  871 1 8-6008 

USAF  ROME  LAB  TECH 
CORRIDOR  W  STE  262  RL  SUL 
26  ELECTR  PKWY  BLD  106 
GRIFFISS  AFB  NY  13441-4514 

AFMC  DOW 

WRIGHT  PATTERSON  AFB  OH  45433-5000 

US  ARMY  FIELD  ARTILLERY  SCHOOL 
ATSF  TSM  TA 
FT  SILL  OK  73503-5600 

US  ARMY  FOREIGN  SCI  TECH  CTR 
CM 

220  7TH  STREET  NE 
CHARLOTTESVILLE  VA  22448-5000 
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NAVAL  SURFACE  WEAPONS  CTR 
CODE  G63 

DAHLGREN  VA  22448-5000 

US  ARMY  OEC 

CSTE  EFS 

PARK  CENTER  IV 

4501  FORD  AVE 

ALEXANDRIA  VA  22302-1458 

US  ARMY  CORPS  OF  ENGRS 
ENGR  TOPOGRAPHICS  LAB 
ETL  GS  LB 

FT  BEL  VOIR  VA  22060 

US  ARMY  TOPO  ENGR  CTR 
CETEC  ZC  1 

FT  BEL  VOIR  VA  22060-5546 

SCI  AND  TECHNOLOGY 
101  RESEARCH  DRIVE 
HAMPTON  VA  23666-1340 

US  ARMY  NUCLEAR  CML  AGCY 
MONA  ZB  BLDG  2073 
SPRINGFIELD  VA  22150-3198 

USATRADOC 
ATCD  FA 

FT  MONROE  VA  2365 1-5170 

ATRC  WSS  R 
WSMRNM  88002-5502 

ARMY  RESEARCH  LABORATORY 
AMSRL  IS  S 
INFO  SCI  &  TECH  DIR 
WSMR  NM  88002-5501 

ARMY  RESEARCH  LABORATORY 
AMSRL  IS  E 
INFO  SCI  &  TECH  DIR 
WSMR  NM  88002-5501 

ARMY  RESEARCH  LABORATORY 
AMSRL  IS  W 
INFO  SCI  &  TECH  DIR 
WSMRNM  88002-5501 

DTIC 

8725  JOHN  J  KINGMAN  RD 
STE  0944 

FT  BEL  VOIR  VA  22060-6218 

US  ARMY  MISSILE  CMND 
AMSMI 

REDSTONE  ARSENAL  AL  35898-5243 


US  ARMY  DUGWAY  PROVING  GRD 
STEDP3 

DUGWAY  UT  84022-5000 

USTRADOC 
ATCD  FA 

FT  MONROE  VA  23651-5170 

WSMR  TECH  LIBRARY  BR 
STEWS  IM  IT 
WSMR  NM  88002 

ARMY  RESEARCH  LABORATORY 
AMSRL  IS  EA 
BATTLEFIELD  ENVIR  DIV 
ATTN  D  BROWN 
WSMR  NM  88002-5501 

ARMY  RESEARCH  LABORATORY 
AMSRL  IS  EA 
BATTLEFIELD  ENVIR  DIV 
ATTN  E  MEASURE 
WSMR  NM  88002-5501 

ARMY  RESEARCH  LABORATORY 
AMSRL  IS  EA 
BATTLEFIELD  ENVIR  DIV 
ATTN  J  COGAN 
WSMR  NM  88002-5501 

Record  copy 


TOTAL 


